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We theoretically study the formation of Brillouin precursors in Debye media. We point out 
that the precursors are visible only at propagation distances such that the impulse response of the 
medium is essentially determined by the frequency dependence of its absorption and is practically 
Gaussian. By simple convolution, we then obtain explicit analytical expressions of the transmitted 
waves generated by reference incident waves, distinguishing precursor and main signal by a simple 
examination of the long-time behavior of the overall signal. These expressions are in good agreement 
with the signals obtained in numerical or real experiments performed on water in the radio-frequency 
domain and explain in particular some observed shapes of the precursor. Results are obtained for 
other remarkable incident waves. In addition, we show quite generally that the shape of the Brillouin 
precursor appearing alone at sufficiently large propagation distance and the law giving its amplitude 
as a function of this distance do not depend on the precise form of the incident wave but only on 
its integral properties. The incidence of a static conductivity of the medium is also examined and 
explicit analytical results are again given in the limit of weak and strong conductivities. 

PACS numbers: 42.50.Md, 41.20.Jb, 42.25.Bs 


I. INTRODUCTION 

In their celebrated papers EH3 on the propagation of a 
step-modulated optical wave through a single-resonance 
Lorentz medium, Sommerfeld and Brillouin found that, 
in suitable conditions, the arrival of the signal at the car¬ 
rier frequency (“main signal”) is preceded by that of a 
first and a second transient, respectively generated from 
the high and low frequencies contained in the spectrum 
of the incident wave. More than a century after their 
discovery, these transients, currently named Sommerfeld 
and Brillouin precursors, continue to raise considerable 
interest, due in particular to their sub-exponential at¬ 
tenuation with the propagation distance. An abundant 
bibliography on the subject is given in Ref [ 51 an d more 
recent related studies are reported in Refs [6Ml0l| . 

In fact the simultaneous observation of well- 
distinguishable Sommerfeld and Brillouin precursors, as 
considered by their discoverers, requires experimental 
conditions that seem unrealizable in optics A quali¬ 
tative demonstration of separated Sommerfeld and Bril¬ 
louin precursors has only been performed in microwaves 
by using guiding structures with dispersion characteris¬ 
tics similar to those of a Lorentz medium El- In op¬ 
tics, only the unique precursor resulting from a complete 
overlapping of the Sommerfeld and Brillouin precursors 
HI . fl3| has been actually observed (TjI i~i~9l | . 

Precursors are obviously not specific to electromag¬ 
netic waves and Lorentz media. A convincing demon¬ 
stration of Sommerfeld precursors has been performed by 
using elastic waves propagating on a liquid surface [20]. 
We consider here electromagnetic waves propagating in 
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dielectric media whose susceptibility is determined by the 
partial orientation of molecular dipoles under the effect 
of the applied electric field [2l|. These currently called 
Debye media are opaque (transparent) at high (low) fre¬ 
quency. Only the Brillouin precursor is thus expected 
to be observable in good conditions, without the over¬ 
lapping problem encountered in Lorentz media. Prop¬ 
agation of waves in Debye media has been extensively 
studied in the past, with particular attention paid to 
the important case of water in the radio-frequency and 
microwave domains. For papers related to the present 
study, see, e.g., Refs [l(| I22| - l34| . Albanese et al. were 
the first to refer to the Brillouin precursors in Debye me¬ 
dia [22]. They numerically studied the propagation in 
water of sine waves at 1 GHz with square or trapezoidal 
envelopes and noticed “a well-formed transient field that 
appears similar to the Brillouin precursor observed in me¬ 
dia with anomalous dispersion”. Stoudt et al. [27) per¬ 
formed the corresponding experiments, achieving direct 
detection of Brillouin precursors. Indirect experimental 
demonstrations were achieved later (30[ [32) . Oughstun 
et al. theoretically studied at length the propagation of 
waves with a step, rectangular or trapezoidal envelope 
in Debye media by combining numerical computations 
and analytical calculations using saddle-point methods 
sm, [341]. Important features of the transmitted wave 
and, in particular, some remarkable shapes of the Bril¬ 
louin precursor were unfortunately overlooked in these 
works. 

In the present paper, the problem is greatly simpli¬ 
fied by remarking that the precursor is only visible (dis¬ 
cernible from the main signal) at propagation distances 
for which the impulse response of the medium is practi¬ 
cally Gaussian. By convolution procedures, this enables 
us to obtain fully analytical expressions of the transmit¬ 
ted waves generated by different incident waves, distin- 
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guishing main signal and precursor by a simple exam¬ 
ination of the long-time behavior of the overall signal. 
The arrangement of our paper is as follows. In Sec El 
we present our approximation and give the correspond¬ 
ing impulse response of the medium. We derive in Sec lIIIl 
the transmitted wave generated from incident sine waves 
with a step or rectangular envelope. The response to 
sine-waves with linearly varying amplitude is studied in 
Sec 030 General properties of the precursors in the strict 
asymptotic limit are established in Sec [V] and the effects 
of a static conductivity of the medium are examined in 
Sec I VII We finally conclude in Sec IVIII by summarizing 
our main results. 


II. IMPULSE RESPONSE OF THE MEDIUM 

As in the experiments on water reported in Refs 
[ril ;27], we consider transverse electromagnetic waves 
propagating in a coaxial transmission line containing the 
Debye medium. This coaxial geometry has the advan¬ 
tage of having a_ “flat frequency response down to and 
including DC” 0, as required for a good observation of 
the Brillouin precursors. We denote £ as the length of 
the transmission line, u(0,f) [v(£,t)] as the voltage at its 
input [output] inside the medium and V (0, w) [ V(£, w)] as 
its Laplace-Fourier transform. In the frequency domain 
the medium is fully characterized by its transfer function 
H(£,uj) relating V(£, w) to V(0, w) [35l] : 


V(£,w) = H{£,lo)V(0,uj) 


with 


H(£, u>) = exp 



(1) 

( 2 ) 


where to, c and n(u>) are the (angular) frequency, the light 
velocity in vacuum and the complex refractive index of 
the medium, respectively. For Debye media, we have 


In the time domain, the medium will be character¬ 
ized by its impulse response h(£,t), which is the inverse 
Fourier transform of H(£,oj), and the output voltage is 
given by the convolution product 


v(£, t) = h(£, t) (g) v(0,t) (4) 

h(£, t) has no exact analytical form. Some general 
properties of h(£, t) can however be derived from 
H(£,lo). The relation f+V h(£,t) dt = H(£, 0) = 1 shows 
that it has a unit area. The location of its center of 
gravity t B = t h (£, t) dt) / h(£, t ) dt ), 

its centered second moment or variance 
cr 2 = - t B ) 2 h(£, t ) dt) / h{£, t.) dt) 

and its centered third moment /i 3 = 

- ts) 3 h{£, t) dt) / (/+“ h(£, t) dt) are re¬ 
spectively equal to the cumulants ki , &2 and fc 3 of 
the transfer function [37j] as defined in the following 
expansion: 


H (£, oj) = exp 



(~^r 

n\ 



(5) 


In the case of a Debye medium, we deduce from Eqs.© 

and ([3j t B = n 0 £/c , a 2 = (^) and /r 3 = 

3 ( n o- n L)( 3n o+ n L) f tr 2 \ , ., , 

^ I These expressions show that 

the center of gravity of h(£, t) propagates at the phase 
velocity at zero frequency (equal to the group velocity 
for this frequency) and has a root-mean-square duration 


a proportional to \/7 and a positive skewness or asym- 

c — v 3 _ 3 ( 3n 0 +"L) rcr 

metry t, - - , Vi' 

4 V n o( Jl o- n ?o) 

The previous results are valid for arbitrary propagation 
distances. The expression of the skewness (£ oc 1 /Vi) 
suggests that the expansion of Eq. flSJ| may be limited to 
the term n = 2 when £ is large enough. Taking the origin 
of time at t = t B , the transfer function and the impulse 
response are reduced to the Gaussians 


h(ui) 



n 0 n oo 

1 + liOT 


(3) 


H(£, u) = = e ~ a( “ )£ (6) 


In this expression n oq, no an d r respectively denote the 
refractive index at high frequency (Debye plateau), the 
refractive index at vanishing frequency and the orienta¬ 
tion relaxation time of the polar molecules. Eq.([3]) pro¬ 
vides a good approximation of the refractive index of De¬ 
bye media in the radiofrequency and microwave domains. 
For deionized water, often taken as reference, ~ 79 , 
71 00 ~ 5.5 and r « 8.5 ps may be considered as typical 
values {27] and are used in our calculations. For the fre¬ 


quencies up to 100 GHz the real and imaginary parts of 
the refractive index derived from Eq. (O with these pa¬ 
rameters fit very well the measurements reported in Ref 


h (£, t ) = JLe - f * 2 (7) 

where a(u>) is the absorption coefficient of the medium 
at the frequency ui and 


P = 


1 

oVi 


cn 0 


2 £t {nl-nlc) 


( 8 ) 


h(£,t) has then a peak amplitude (a duration) propor¬ 
tional to l/y/I (VI) with a unit area, as expected. It 
meets the principle of relativistic causality [38| as long 

as exp —/3 2 (n 0 — noo) 2 ( £/c ) 2 is negligible, a condition 
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superabundantly satisfied for the propagation distances 
at which the Brillouin precursor is visible. Decompos¬ 
ing the medium in m subsections of impulse response 
h(£/m,t), the Gaussian form of h(£,t) given by Eq. 0 
may be considered as a consequence of the central limit 
theorem in a deterministic case [35]. It is also the limit 
when (Hb 1 of the expressions of h(£, t) obtained by 
direct studies in the time-domain [23|, (25j, [28|. When 
£ <C 1, the condition fits 1 is automatically fulfilled 
and the Gaussian form h(£,t) given by Eq. 0 is ex¬ 
pected to hold. This is illustrated Fig [1] obtained for 
£ = 10 cm in deionized water. We have then (3tB ~ 13.7 



Figure 1: (Color online) Impulse response of deionized water 
as a function of time for a propagation distance l = 10 cm 
leading to ts = 2.96 ns and /3 = 4.63 ns -1 . The solid [dashed] 
line is the exact numerical result obtained by FFT [the Gaus¬ 
sian shape given by Eq. 01 . Inset: the same for £ = 1 mm, 
leading to ts = 29.6 ps and /3 = 46.3 ns -1 . 

and the result given by Eq. 0 is actually very close 
to the exact result, numerically derived from the exact 
transfer function by fast Fourier transform (FFT). The 
main effect of the residual skewness is a very slight re¬ 
duction of the time of the maximum by about 9.5 ps 
whereas ts ~ 2.96ns. For £ = 1 cm (fits ~ 4.3), the 
skewness increases and the rise of h(£,t) becomes signif¬ 
icantly steeper than its fall. However the Gaussian form 
of Eq. 0 remains a reasonable approximation of the 
exact result. This approximation completely fails when 
pts 1 ■ The impulse response then begins at the time 
n 00 £/c {38 ] by a Dirac peak and a discontinuity. The 
weight w of the former and the amplitude Ah of the lat¬ 
ter are easily determined from an asymptotic expansion 
of H(£,uj). We get w = exp [— (rig — n^) £/ (2nooCr)] 
and Ah = w£ (rig — n^) (rig + 3n;( 0 ) / (8CT 2 nj^). The 
inset of Fig. [T] shows the impulse response obtained for 
£ = 1 mm {fits ~ 1.37). We have then w ~ 2.14 x 10” 3 
and Ah « 6.7 ns” 1 . The previous results are only given 
for completeness. Indeed, it will be shown later that 
the Brillouin precursor emerges from the main signal for 


£ > 20 cm (0tB > 19.4 ) and the Gaussian approxima¬ 
tion of h(£, t) is then excellent. Note that this remarkable 
form of the impulse response is not specific to the Debye 
medium but holds at large enough propagation distances 
whenever the transfer function can be expanded in cumu- 
lants. It applies in particular to polar media with relax¬ 
ation mechanisms more complex than those considered 
here. We incidentally note that the addition of a second 
relaxation time in the Debye model as considered in Refs 
sm does not significantly modify the parameters ts 
and f3. 

The considerable advantage of the Gaussian impulse 
response is that it can be convoluted with a great num¬ 
ber of input signals to provide exact analytical forms 
of the output signal. This is the case, e.g., when 
v(0, t) is a step, an algebraic function, an exponential, 
a Gaussian or an error function, when these functions 
modulate a sine wave (eventually linearly chirped) or 
when u(0, t) is a combination of the previous signals 
d 0, [H IH EE Hi Hi, HU- The relation V(£, 0) = 
H(£, 0)R(0, 0) quite generally shows that the area of the 
output signal will be always equal to that of the input sig¬ 
nal. We also remark that the Gaussian impulse response 
is obtained by neglecting the effects of the group-velocity 
dispersion. This explains why the shape of the output 
signals observed in the experiments reported in Ref [27j 
was well reproduced by numerical simulations taking only 
into account the frequency-dependence of the medium 
absorption. This also means that the Brillouin precur¬ 
sors actually observed or observable in Debye media es¬ 
sentially originate from the latter and that group-velocity 
dispersion plays no significant role in their formation. 

In the previous theoretical analysis it is assumed as 
usual that the input and output signals are measured 
inside the Debye medium. This is generally not the case 
in the experiments. H(£,uj) should then be multiplied by 
the transfer function T(u>) taking into account the losses 
at the air-medium and medium-air interfaces, that is 


2 

X 

2 

1 + n(w)_ 

1 + l/n(w)_ 



4n(w) 

[i + «M ] 2 


(9) 


The main effect of T(u>) will be an overall reduction of 
the amplitude of the output voltage by a factor 1/T(0) = 
(no + l) 2 / (4ng) that is 2.75 for deionized water. The 
frequency-dependent effects can be determined by ex¬ 
panding T(w)/T( 0) in cumulants as made for H(£,uj) 
and exploiting the additivity property of the cumulants. 
Again for deionized water, we find that ts (cr 2 ) is re¬ 
duced by about 3.15ps (15.7ps 2 ). For £ > 10cm, these 
quantities are fully negligible compared to those associ¬ 
ated with H(£,lo). The output signals determined from 
H(£,ui) in the following should thus be simply multiplied 
by T(0) = 4no/ (1 + no) and delayed by the extra tran¬ 
sit times in air when the voltages are measured outside 
the medium. 
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III. RESPONSE TO STEP-MODULATED 
SINE-WAVES 

We consider in this section input signals of the form 
unit) sin (oj c t) or unit) cos (u) c t), where unit) is the 
Heaviside unit step function. They consist of a sine 
wave of angular frequency lo c (the “carrier”) switched 
on at time t = 0. The former is the “canonical” signal 
used by Sommerfeld and Brillouin. Both cases can be 
jointly studied by considering the complex input signal 
v(0, t) = un(t) exp (iuj c t). Its convolution with the Gaus¬ 
sian impulse response gives for the corresponding output 
signal 


v(£, t) « - [1 + erf (fit + iy)\ e v e zulct (10) 

Here erf(xr) indicates the error function, t is the time 
retarded by tg, and y = uj c / (2/3) = \/a c £, where a c is a 
short-hand notation for a (w c ). Taking ts (1//?) as time 
origin (time scale), Ea. TOH) shows that the output signal 
only depends on y. For sufficiently large t, this signal 
becomes 


v m {l, t) » i [1 + erf (fit)] (11) 

and tends to v rn (£, t) « e -a c l e %u c t w j ien ^ —> oo. This 
part of v(£,t) given by Ea. dlll) may naturally be identi¬ 
fied to the main signal. The Brillouin precursor is then 
given by the remaining part, that is 


(£, t) ~ — erf ^/3 1 + iyja c £^ —erf (fit) 


—ot c lAoj c t 


( 12 ) 

Since ui c = 2/3\Ja c £ , all the previous signals appear as 
universal functions of fit or of uj c t that only depend on the 
optical thickness a c £ of the medium at the frequency uj c 
of the carrier. To be definite we consider in the following 
deionized water and a carrier frequency u c = 27rx 10 9 s -1 , 
i.e., a period Th = Ins as often considered in the liter¬ 
ature [5f, 122, [271, [29|. We, however, emphasize that our 
analytical results are quite general. 

The responses to the canonical input signal 
un(t) sin (u> c t) are obtained by taking the imaginary part 
of Eqs. m)-m- Figure [3] shows the output signals 
v(£,t) obtained for different values of a c £. As expected, 
they reproduce very well the exact signals obtained by 
a FFT calculation of the inverse Fourier transform of 
H(£,(jj)V(0,uj). For a c £ =!(•£« 21.6cm), the precursor 
appears as a small overshot on v(£,t), which is hardly 
visible for smaller optical thicknesses. For a c £ = 3, we 
have an example of well developed precursor dominating 
the main signal. Finally for a c £ = 10, the main signal is 
extremely small and only the precursor remains visible. 
Figure [3] shows the corresponding precursors alone, as 
defined by Eq. m- Owing to the symmetry properties 
of the real and imaginary parts of the error function, 
they are even functions of t. Their peak amplitude is 



Retarded time (ns) 

Figure 2: (Color online) Responses to the input signal 
uh (t) sin(uj c t) as a function of the time retarded by ts, for 
a c £ = 1, 3 and 10 (from top to bottom). The solid (dashed) 
lines give the exact numerical solution (our analytical solu¬ 
tion in terms of error functions). Note that for the carrier 
frequency considered (w c /(27r) = 1 ns -1 ), the optical thick¬ 
ness a c £ = 1 is attained for i = 21.6 cm. 


simply 


A p (£) = Im [v p (£, 0)] = crf (13) 

whereas their area S p (£) , also deduced from Eq. (fl^l) , 
reads as 


S P (£) 


(1 _ e -2 
UJ C 


(14) 


Remarkably enough, their shape is practically Gaussian. 
The precursor part of the output signal can thus be writ¬ 
ten under the approximate form 

v p (£, t) = A p (£) e-PZ t2 , (15) 


where the parameter f} p is easily determined by combin¬ 
ing Eqs. m and (1141) . We get 


Pp = 


yj 7T a c £ erf (iy/a c £) 

i (g a-ct - g—a c t) 


P- 


(16) 


As shown in Fig. [3l Eq. (fl5l) . with A p and j3 p given 
by Eq. m and (USD, provides an excellent approxima¬ 
tion of the exact result. It strictly gives the exact result 
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Retarded time (ns) 

Figure 3: (Color online) Analytical forms of the precursor 
contribution to the signals of Fig. [2] (solid lines) and their ap¬ 
proximation by the Gaussian given in Eq. (1151) (dashed lines). 
Inset : ratios of ft p (solid line) and A v (dashed line) over their 
asymptotic values as a function of the optical thickness. 



Retarded time (ns) 

Figure 4: (Color online) Same as Fig. [2] when a(0, t) ■ 
un{t) cos (uj c t) for a c £ = 3, 5, and 10 (from top to bottom). 


when a c £ —> oo. From the asymptotic form of the error 
function 1391. we then get A v (£) = - 1 , = —and 

1 - ,J ’ ° / 2y/na c i UcV 1 

ftp = ft, that is, 


Vp {£, t ) 



e -/3 2 t 2 

2^/na c £ 


(17) 


The well-known law A p {£) oc 1 /\f£ is retrieved but we 
remark that it only holds for very large optical thick¬ 
nesses. When a c £ is only large (such that e~ act <C 1), it 
is possible to obtain a better approximation of v p (£, t) by 
considering one more term in the asymptotic expansion 
of the error function [39]. This leads to a simple multipli¬ 
cation of ft by y/l + 1/ (a c £) in the first form of Eq. m- 
The inset of Fig. [3] more generally shows how the ratios 
of ftp and A p over their asymptotic limit (respectively ft 
and ^ ^p ) vary as a function of a c £. 

We examine briefly the case where i>(0,f) = 

un(t) cos (u> c t). Such an input field originates a class of 
precursors whose amplitude depends on £ according to 
a law differing from the previous one. We incidentally 
remark that, for uj c = 0 (no carrier), the output signal is 
nothing other than the step response of the medium [35| 
that reads as a(£, t) = [1 + erf {ft f)] /2, a result in agree¬ 
ment with the observations reported in Refs Hill. For 
oj c ^ 0 , the output signals are obtained by taking the 
real parts of Eqs. (fII)l) - (fT^l) . Figure [I] shows the output 
signals v(£,t) obtained for a c £ = 3, 5, and 10. The pre¬ 
cursor appears for propagation distances larger than in 
the previous case and is not discernible from the main 
signal when a c £ = 1. The output signals v(£,t) derived 


from Eq. (HUD again fit very well the exact signals ob¬ 
tained by FFT. The contribution v p {£, t) of the precursor 
to the output signal (Fig. [5j) is now an odd function of 
t. Its slope at t = 0 can be determined by expanding Eq. 
(fl^l) at the first order in t. We get v p {£, 0) = — B p , with 


B p = ft 


z-erf (iy/a c £) „ 

- V v j e _ ac£ 


1 — e 


>0. 


(18) 

The precursors are now well fitted by derivatives of Gaus¬ 


sian (Fig. [5]) and can be written as v p {£, t) « —B p te~7p t . 
In fact, 7 p , achieving the best fit, never considerably dif¬ 
fers from ft p . The precursor takes then the approximate 
form 


v p (£,t) « —B p te PA 


B p d — Ppt 2 

2 ft% dt 6 


with a peak amplitude 


A p {£) 


B p 

PpV2e 


(19) 


( 20 ) 


As shown in Fig. [5J Eq. (fTTil) provides a satisfac¬ 
tory approximation of the precursor. It gives the ex¬ 
act result when a c £ —> oo. In this limit, ft p = ft, 
B p = 2ft 3 / (ui/y/TT^j, A p {£) = 1/ {2a c £y/2eir) and 


v p (£,t) 
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Retarded time (ns) 

Figure 5: (Color online) Analytical forms of the precursor 
contribution to the signals of Fig. U (solid lines) and their 
approximation by the Gaussian derivative given Eq. m 
(dashed lines). Inset is the same that of Fig. [3] 

Since w c ii#(f) cos (u c t) = jp [unit) sin (w c t)] and 
v p (£,t) = v(£,t) when a c l oo, this precursor is the 
time derivative of that obtained in the canonical case 
[see Ea. (fT71) | divided by uj c . Its amplitude now scales as 
1/i instead as 1/in the previous case. As shown in the 
inset of Fig. 0 this law only holds for very large optical 
thicknesses. 

From the results obtained with the input signals 
unit) sin (ui c t) and un(t) cos (ui c t), the responses to the 
more general input signal unit) sin (w c f + tp) are easily 
determined. The precursor is then a superposition of a 
Gaussian and a Gaussian derivative with relative ampli¬ 
tudes depending on tp and on the propagation distance. 

In real or numerical experiments, the input signal 
has obviously a finite duration. One generally consid¬ 
ers input signals with a rectangular envelope of dura¬ 
tion T long enough to avoid having the precursors gen¬ 
erated by the rise and the fall of the pulse overlap. 
The corresponding complex signal reads as vn(0,t) = 
[u H (t) — unit — T)[ e lulct , that is 

v fl (0, t) = u H (t ) e iUct - e iuicT u H (t - T ) e iu ^ t ~ T '> . (22) 

It generates the output signal 

vr(£, t) = v(£, t ) — e*“° T v(£, t — T), (23) 

where v(£,t) is given by Eq. (flCil) . No matter T, the 
mains signals generated by the rise and the fall of the 
pulse destructively interferes when 0 (t — T) \fa c £- 
On the other hand, the second precursor (a postcursor 
since it follows the main signal) generally differs from the 
first one owing to the phase factor e lu>cT . Precursor and 
postcursor are identical only when T = (2 n + 1) T c / 2, 
where n is an integer. They have the same shape and 


amplitude but opposite signs when T = nT c . Both 
cases have been evidenced in the experiments reported 
in Ref na. The largest difference between precursor 
and postcursor is attained when T = (2 n+ 1) T c /4. In 
this case e lLJcT = (—1 ) n i. For an input signal v(0,t) = 
[unit) — un{t — T)] sin (ui c t), we get an output signal 

vr(£, t ) = Im [v{£, t)] - (-1)" Re [v{£, t - T)] (24) 

The precursor is then Gaussian [see Eq. (fT5l) ] whereas 
the postcursor is a Gaussian derivative [see Eq. dl9ll ] 
of smaller amplitude. The opposite occurs for an input 
signal u(0, t) = [un(t) — un{t — T)] cos (u c t). These var¬ 
ious behaviors are illustrated in Fig. [Gj 



Figure 6: (Color online) Responses to input signals with a 
long rectangular envelope for a c £ = 5. The upper curves are 
obtained for v(0,t) oc sin(u;ct) with T = 10 If (solid line) 
and T = 10.5 T c (dashed line). The lower curves are ob¬ 
tained for T = 10.25 T c with v(0, t ) oc sin(ai c t) (solid line) 
and u(0, t) oc cos (u] c t) (dotted line). Note that, in the latter 
case, the postcursor can completely dominate the precursor. 

Equation (|23D remains valid when the two precur¬ 
sors overlap and interfere to give a unique transient 
that may be considered as a generalized precursor. A 
remarkable behavior is obtained for an input signal 
[unit) — Uff(t — T c )[ cos (w c t). A good overlapping of the 
precursors is achieved when 1/0 T c , that is in the 
asymptotic limit. We then get 

v(£,t) « v p (£,t)-v p (£,t-T c ) « T c v p ^£,t- , (25) 
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where v p (£,t) is given by Eq. (lUTl) . and finally 


( 2/3v^ 

1 — (e~P 202 ) 

1 “c J 

d6' 2 V ) 


- A oo (l-2p 2 e 2) ^e-^ 0 \ (26) 

where 0 = t — T c /2 and = (yUr/2) (a c £)~ 3 ^ 2 . The 
output signal has now the shape of a Gaussian second- 
derivative centered at t = T c /2 and a peak amplitude Aoo 
scaling as t' -3 / 2 . For any optical thickness, it is easily 
shown from Ea. (l23l) that v(£, 6) is an even function of 0 
with a minimum at 8 = 0 of algebraic amplitude —A, 
with 



A = e 


that tends to Aoa 
Eq. (IU51) that 


Re 


erf 


f 

V W c 


when a c £ —> oo. 



(27) 


It also results from 


Figure 7: (Color online) Response to the input signal 
[uH(t) — uh (t — T c )] cos(wct) for a c £ = 5. The solid, dashed 
and dotted lines are the exact FFT solution, the approximate 
analytical solution given Eq. ([29jl. and the asymptotic solu¬ 
tion given Eq. (12611 . respectively. Inset: input signal. 


v(e,9) 



erf 


( 27T/3 . UJ C \ 

\— + 'w) 


(28) 


Proceeding as in the previous section, we get for the main 
signal and the precursor 


for 9 = ±7r/w c , that is, at the switching times of the input 
signal delayed by ts- Unexpectedly enough, extensive 
numerical simulations show that v (£, 9) is very well fitted 
by a Gaussian second derivative as soon as a c £ > 1. The 
output signal then reads as 


_ q ^ 

= -yyjr {{pt + [1 +erf (/?i)] j e lulc \ 

(31) 

v p (£,t) = v(£,t) - v m (£,t). (32) 


v (£, 6) = - A (1 - 2t f/3 2 9 2 ) e-^P 262 , (29) 

where A is given by Eqs. (E71) and r] is a parameter ob¬ 
tained by combing Eq. (l27l) and (128)1 . Figure [7] shows 
the result obtained by this method for a c £ = 5. It per¬ 
fectly fits the exact numerical result obtained by FFT. 
The asymptotic result of Eq. (1^1) is also given for com¬ 
parison. In spite of the moderate value of the optical 
thickness, it is reasonable approximation of the exact re¬ 
sult. 


IV. RESPONSE TO SINE-WAVES WITH 
LINEARLY VARYING AMPLITUDE 


We consider in this section input signals of the form 
( t/T r ) unit) sin (w c f) or ( t/T r ) unit) cos (u> c t) where T r is 
the time for which their amplitude equals 1 (rise time). 
Such signals can generate precursors preceding signifi¬ 
cantly the rise of the main signal. The corresponding 
complex signal reads as v(0,t) = (t/T r ) UH(t) e zulct . By 
convoluting this signal with the Gaussian impulse re¬ 
sponse, we get 




e -0 2 t 2 e -a c e 
2/3T rV / 7r + 2 /3T r X 




1 + erf 




]} 


e iu c t 


(30) 


All these signals are inversely proportional to T r , and 
w c being fixed, it is convenient to normalize them to 
1 / (u c T r ). 

Figure [8] shows the output signal v(£,t) = lm[v(£,t)] 
generated by v(0,t) = (t/T r ) un{t) sin (oj c t) for different 
optical thicknesses. The precursor becomes discernible 
from the main field for a c £ > 6 and is fully separated from 
it for a c £ = 12. In every case the analytical result derived 
from Eq. (1301) fits very well the exact numerical result 
obtained by FFT. The contribution v p (£, t) = Im [v p {£, £)] 
of the precursor to the output signal is now an odd func¬ 
tion of the retarded time t and, for a c £ > 6 , is well fitted 
by a Gaussian derivative (Fig. [9]). This result is rigorous 
when a c £ —> oo . By means of asymptotic expansions 
of the error function limited to its first non zero term in 
1 /y/aj., we then get 



the precursor amplitude scaling as 1 /£. 

When u(0, t) = ( t/T r ) un{t) cos (w c f) , the overall out¬ 
put signal and the precursor read as v(£,t) = Re [v(£,t)] 
and v p (£,t) = Re[u p (f', £)] . The precursor becomes dis¬ 
cernible from the main field for a c £ > 4.5 (Fig. [TUI) and, 
again owing to the symmetry properties of the error func¬ 
tion, is an even function of t. At t = 0, v p (£,t) has the 























Retarded time (ns) 



Retarded time (ns) 

Figure 10: (Color online) Same as Fig. [8] with v(0,t) = 
(t/T r )uH (t) cos(cj c t) for a c (. = 4.5, 7, and 12 (from top to 
bottom). 


get from Eq. (1551) 


Figure 8: (Color online) Responses to the input signal 
(t/T r )uH(t) sin(cj c t) normalized to l/(a i c T r ) for a c £ = 6, 8, 
and 12 (from top to bottom). The solid (dashed) line is the 
FFT exact solution (the approximate analytical solution in 
terms of error functions). 



Retarded time (ns) 

Figure 9: (Color online) Analytical forms of the precursor 
contribution to the signals of Fig. [8] (solid lines) and their 
best fit by Gaussian derivatives (dashed lines). 


remarkable value 


_ aJ e rftn/33^^ 

i 

(34) 

where A > 0 as soon as a c i > 1. For a c £ > 4.5, the 
precursor is well fitted by a (negative) Gaussian (Fig. 
EH). This solution is exact when a c i —> oo and we then 



0 (£,t) = - 


P 




= -A n 


/7T 


(35) 


where A^ = 1/ (2ui c T r y/Tra c t\ is the limit of A when 
a c l —> oo and scales as 1 j\fl. 



Retarded time (ns) 

Figure 11: (Color online) Analytical forms of the precursor 
contribution to the signals of Fig. ESI (solid lines) and their 
best fit by Gaussians (dashed lines). 

When the precursors are generated by a single discon¬ 
tinuity of the input signal, their shape in the asymp¬ 
totic limit strongly depends on the order of this dis¬ 
continuity. Equations (E71) . (f2T1) . (1551) . and (l35l) lead us 
to conjecture that the precursor has a Gaussian (Gaus¬ 
sian derivative) shape when the discontinuity order is 
odd (even). Complementary calculations made when 
v(0, t) = (t/T r ) n unit) sin (w c f) or ( t/T r ) n unit) cos (w c f) 
for n = 2 and n = 3 support this conjecture but it should 
be remarked that the precursor amplitude is a rapidly de¬ 
creasing function of n . 
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The effect of the rise time of the input signal on the 
precursors is generally studied by considering input sig¬ 
nals whose amplitude linearly increases du ring a time 
T r < T c and is then maintained constant !'i I22l.l27l.i34l]. 
The corresponding complex signal may be written as 


VL{0,t) = ~~u H (t) e^ ct 
1 r 

_ e iu lc Tr t _zJL UH ^- T r)e iu ^t-Tr)_ (36 ) 
1 r 

It generates the output signal 

v L (t,t) = v(£,t) - e iu ^v{£,t- T r ), (37) 

where v(£, t ) is given by Eq. (1501) . For j3(t — T r ) a c £, 
Eq. (1571) is reduced to vl(£, t) = e _ “ c€ e* Wct , as expected 
since the long-term behavior of the main signal should 
not depend on the rise time of the input signal. Remark¬ 
able behaviors are obtained in the asymptotic limit when 
T r = T c and T r = TJ 2. When T r = T c , the precursors 
generated by the first and second discontinuities of the 
slope of the input signal interfere nearly destructively 
and the output signal are again given by Eq. (l25ll where 
v p (£,t) is now given by Eq. (1551) for v(0,t) oc sin (oj c t) 
and by Eq. (1551) for v(0,t) oc cos (u c t) . The output 
signals then reads as 




2/3 d 2 


(e-^) 


_ (1 - 2 (3 2 6 2 ) e-P 2 
2 ( a c £) 3 ^ 2 

(amplitude oc £ -3 / 2 ) in the former case and as 

d ( p6e-P 2(>2 


v(t,0) = - 


d6 


(e-? 2 * 2 ) 


2a c £^ 


(38) 


(39) 


(amplitude oc l/£) in the latter case. In both cases, 
6 = t — T c /2 . When T r = T c /2 , the two precur¬ 
sors interfere constructively. For u(0 ,t) oc sin (ui c t) and 
u(0, t) oc cos (w c t), we respectively get 


v(i,e) = 


4/3 


w 2 7r 3 / 2 dd 
(amplitude oc l/£) and 


(e^ 202 ) 


2/3 6e~ 




W C 7T' 


3/2 






7r 3 / 2 ^/ a c l 


(40) 


(41) 


(amplitude oc ^ _1 / 2 ), with 6 = t — T c /4 in both cases. 
Figure [T5] shows the output signals given by Eqs. (1571) 
combined with Eq. (1501) in the four cases considered in 
this paragraph. The optical thickness aJL = 7 has been 
chosen in order that the main signal is dominated by 
the precursor but remains visible. The precursors ob¬ 
tained at the asymptotic limit are given for reference. 


Though the asymptotic limit is not attained, they pro¬ 
vide a reasonable approximation of the complete signals. 
Complementary simulations obviously show that as a c £ 
is larger, the fit of the output signals by the asymptotic 
form of the precursors is better. Signals as those shown 
Fig. [12] have been actually observed in the numerical ex¬ 
periments reported in Refs @0130 but the fact 
that the precursor shape may be a Gaussian, a Gaus¬ 
sian first derivative or a Gaussian second derivative is 
either not clearly stated or completely overlooked. All 
these numerical simulations were made by using a trape¬ 
zoidal modulation with the same rise and fall times and 
a plateau whose duration is much larger than that of 
the precursors. As in the case of a rectangular modula- 


0.01 

0 

- 0.01 


- 0 


- - 0.05 


-2 0 2 4 -2 0 2 4 

Retarded time (ns) 

Figure 12: (Color online) Responses (solid line) to input sig¬ 
nals with trapezoidal envelope of rise times T r = T c or T c /2 
for a c £ = 7. (a) T r = T c , v(0,t) oc sin(w c t); (b) T r = T c , 
v(0,t) oc cos(w c t); (c) T r = T c / 2, v(0,t) oc sin(u;ct); and (d) 
T r = T c / 2, u(0, t) oc cos (uj c t) . The asymptotic solutions 
(dashed lines) are given for comparison. 



tion, depending on the plateau duration, the precursor 
and the postcursor may have the same or opposite signs 
t22h I23, [34j and even different shapes, e.g. a Gaussian 
associated with a Gaussian first-derivative or a Gaussian 
first derivative associated to a Gaussian second deriva¬ 
tive. We incidentally remark that, by using a modula¬ 
tion with rise time, plateau duration and fall time all 
equal to T c , it would theoretically be possible to gener¬ 
ate a precursor with a Gaussian third derivative shape 
and an amplitude scaling as l/£ 2 . Calculations made in 
this case show that the corresponding precursor would 
emerge from the main field at distances where it would 
be too small to be detected in a real experiment. 
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V. GENERAL ASYMPTOTIC PROPERTIES OF 
THE PRECURSORS 

The shapes of the precursors obtained at the asymp¬ 
totic limit and the dependence of their amplitude with 
the propagation distance (as or £~ 3 / 2 ) are 

not specific to the input signals considered in the pre¬ 
vious sections. As shown below they only result from 
the fact that, in the asymptotic limit, the width of the 
transfer function is infinitely small compared to 

that of the Laplace-Fourier transform V (0, to) of the in¬ 
put signal. The latter can then be developed in cu- 
mulants by keeping only the first of them. It reads as 
V(0,w) = V(0, 0)e _I “ to = Soe~ lcot ° where So is the area 
of the input signal and to its center of gravity [40]. We 
then get V(£, w) = S 0 H{£, Lu)e~ iulto and finally 

v(£,t) = S 0 h(t-t 0 ) (42) 

This result, leading to a precursor amplitude scaling as 
£~ 1 / 2 , is general and applies in particular to the pre¬ 
cursors given by Eqs. (H3, ED, and (El), for which 
the input signals are such that [So, to] = [l/u; c ,0], 
[—1/ (w 2 T r ) ,0] and [—2/ ( 7 rw c ) , T c /4], respectively. 

The previous demonstration obviously fails when the 
input signal has a zero area. We consider in this case the 
first antiderivative of the input signal, namely (Ri;(0, t) = 
f_ v(0, Mi) du 1 , the Fourier transform of which reads as 
F(0,w)/ (iui). Repeating the previous procedure, we get 
V(£,uj) = (iw) uj)e~ lultl and finally 

v(£, t) = S\ h(t — ti) (43) 

where Si and t\ are the area and the center of grav¬ 
ity of Wm(0 ,t) [40]. This general result, leading to a 
precursor amplitude scaling as l/£, applies in particu¬ 
lar to the precursors given by Eqs. ED, HMD, HMD, and 
(HOD- The corresponding input signals are actually such 
that So = 0 whereas [Si, ti] = [l/w 2 ,0], [2/ (u^T r ) , 0], 
[—1/w 2 , T c /2] and [4/ (itlu 2 ) ,T c /4], respectively. 

When So and Si are both zero, the same method 
is applied to the second antiderivative ^v(0,t) = 

f_ S-l 3 m(0, mi) duidu 2 , the Fourier transform of which 

reads as U(0,w)/ (iuj) 2 . Denoting S 2 as the area of 
( 2 )m(0 ,t) and t 2 as its center of gravity (40|, we hnally 
get the output signal 

v{£, t) = S 2 h(t — t 2 ) (44) 

the amplitude of which scales as f? -3 / 2 . This result ap¬ 
plies in particular to the precursors given by Eqs. ED 
and (1381) . The second antiderivative of the correspond¬ 
ing input signals is such that t 2 = T c / 2 in both cases, 
S 2 = 27r/w 3 in the first case and S 2 = 2/w 3 in the second 
case. According to Eq. (l44l) . identical precursors will 
be thus obtained by amplifying the input signal in the 
second case by a factor 7 r. This result is illustrated Fig. 
m which spectacularly shows that quite different input 




Figure 13: (Color online) Responses to the input signals (a) 
(solid line) and (b) (dashed line) for a c i = 12. The asymp¬ 
totic solution of Eq. ED is given for comparison (dotted 
line). 


signals can generate the same precursor as long as they 
have the same integral properties. Figure [13] has been 
drawn for a finite optical thickness ( a c £ = 12) for which 
the precursors have sufficient amplitude and can be really 
observed. Both are close to the asymptotic form given 
by Eq. (1441) from which they become undistinguishable 
for a c £ > 30. 


VI. EFFECTS OF THE MEDIUM 
CONDUCTION 


Debye media have always some conductivity. In the 
reference case of water, this conductivity a is essen¬ 
tially static in the radiofrequency and microwave domains 
[4lL 142] and ranges from about 5 x 10~ 6 S/m for ultra- 
pure water (only due to H + and OH~ ions) to 5 S/m 
for typical seawater (mainly due to dissolved salts). In 
the presence of conductivity, the complex refractive index 
reads as 


n(u>) 



np ~ nlc 

1 + IUJT 


a 

iui£o 


(45) 


When the conductivity is small enough, the main con¬ 
tribution to h(ui) and thus to H(£,uj) is expected to orig¬ 
inate from frequencies such that the last term in Eq. (Uhl) 
may be treated as a perturbation. Repeating the proce¬ 
dure that led to Eq. ©, we then get 


In [H (£, w)] ~ —icots — 


oj 

w 


criZn 


(46) 


2 
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Retarded time (ns) 

Figure 14: (Color online) Impulse response of water for i = 
lm (ts ~ 29.6 ns) and a « 0.1 S/m. The solid (dashed) 
line is the exact numerical result (our approximate analytical 
result). Inset: corresponding response to the canonical input 
signal Uir(t) sin (u> c t) with cj c = 27r x 10 9 s -1 . 


where Zq = \/p,o/ ( n o £ o) is the characteristic impedance 
of the medium at u> = 0 when a = 0. Using again a time 
retarded by ts, we finally get the approximate transfer 
function H a {£,u) 

H a {l,u) = exp e~^ (47) 


previously. As expected, v(£, t ) is fairly well reproduced 
by the analytical expression Im[w(.Rt)] /F, where v(£,t) 
is given by Eq. CUD. and the long tail of h (£, t ) originates 
the upshift of the main signal. 

Before considering the general case, it is instructive 
to examine the effects of the sole conductivity. This is 
achieved by neglecting the second term under the square 
root in Eq. (|45[) . The transfer function is then reduced 
to H c (£,uj) of the form 


H c (£,u) = exp 



(49) 


where s = ioo, too = n-oo£/c and 7 = < 7 / {2eon^). Us¬ 
ing the inverse Laplace transform of exp — too\/s 2 ~ 7 2 

as given in Ref [43] and translating s by 7 , we get the 
impulse response 


h c (£, t) = 7 t c 


u (7 \A 2 -*i>) 




e 7 t u H (t-too) 


+ e lta °5 (t - too) • (50) 


Here R (z) designates the first-order modified Bessel func¬ 
tion and t is the real time. A more rigorous demon¬ 
stration of this result, showing its consistency with the 
boundary conditions, can be found in Ref (44[. For 
t 2 ^ too an d it ^$> 1, Eq. (l50l) takes the asymptotic 
form 


and the approximate impulse response h a (£,t) 

K (^) = ^ e_/32 * 2exp 

At this order of approximation, the effect of the medium 
conductivity is simply to reduce the impulse response 
and thus the response to any input signal by the factor 
F = exp ( at o ° ). The effect of the medium conductivity 
will be negligible when the product of its conductance 
<t£ by its characteristic impedance Zq is small compared 
to 1. For l = lm, this condition is practically realized 
with ordinary drinking water whose typical conductiv¬ 
ity (10 _ 3 S/m) leads to F « 1.02. The value F = 2 is 
attained for a ~ 0.033 S/m and numerical simulations 
shows that h a (£, t) is then a very good approximation of 
the exact impulse response h(£,t). We, however, note 
that the area of h a (£, t ) is only 1/F whereas that of the 
exact impulse function remains equal to 1. The reason 
for this apparent discrepancy is that the low frequencies, 
not fully taken into account in the previous calculation, 
originate a very small but very long tail whose area is 
the missing area. This tail is well visible on Fig. fl4l 
obtained for a « 0.1 S/m. Though F is now as large as 
8.3, Eq. (1481) continue to provide a good approximation 
of the main part of the impulse response. The inset of 
Fig. [ID shows the output signal generated by the canoni¬ 
cal input signal unit) sin (w c f) with ui c = 2n x 10 9 s _1 as 


h c (£, t) 


Pocr £ 2 

47rt 3 


exp - 


/i 0 oT 2 


At 


(51) 


which is maximum for t = t m = no<r £ 2 /6 with an ampli¬ 
tude 


km ~ \[wei {p 0 a£ 2 ) ■ (52) 

Coming back to the general problem, Eq. (H51) is ex¬ 
pected to provide a good approximation of the precursor 
if, at the time of its maximum, the ratio R of h a (£, t) 
over h c (£, t )is large. Using Eqs. (H51l and (I5T1) and tak¬ 
ing into account that Uq n^, we get the approximate 
expression 



In the conditions of the inset of Fig. [H] (£ = 1 m , 
<j w 0.1 S/m.), R was about 14 and Eq. (14811 actually pro¬ 
vided a good approximation of the medium response to 
the canonical input signal with ui c = 2 ttx 10 9 s -1 . Figure 
m shows the result obtained for a « 0.3 S/m, leading to 
R r: 1. The response obtained by considering the sole 
effects of conductivity (4i| approximates fairly well the 
mean value of the exact response for t > ts- The coinci¬ 
dence becomes exact for t ts, the long-time behavior 
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Figure 15: (Color online) Response to the canonical input 
signal for a « 0.3 S/m (other parameters as for the inset of 
Fig. 11411 . The solid (dashed) line is the exact numerical re¬ 
sult (the analytical result obtained by considering the sole 
effects of conductivity). Inset: enlargement of the output sig¬ 
nal around t = ts- Here the dashed line is the analytical 
result derived from Eq. (1481) . 



Time (ns) 

Figure 16: (Color online) Same as Fig. [15] for a « 0.5 S/m. 
Note that the scales of the inset and the main figure are very 
different. 


being mainly determined by the low frequencies where 
the conductivity term prevails on the polarization term 
in Eq. (]45D . On the other hand, the inset shows that the 
beginning of the response is perfectly reproduced by the 
analytical response derived from Eq. (|48p . Main signal 
and precursor obviously decrease with the conductivity. 
For a k, 0.4 S/m, the former remains visible whereas the 
latter appears as a slight overshot of amplitude about 10 5 
times smaller than that of the input signal and would 
be probably undetectable in a real experiment. Figure 
m shows that main signal and precursor become prac¬ 
tically invisible for a ss 0.5 S/m. The overall response 
is well fitted by the analytical function derived from Eq. 


® 0 but Eq. (1551) continue to perfectly reproduce 
the very first beginning of the signal (inset). Finally, as 
the conductivity increases, the fit of the output signal 
by h c (£, t) /u! c improvies. Any trace of precursor then 
disappears and the response is reduced to a broad signal 
of amplitude (duration) scaling as l/£ 2 (as £ 2 ). See Eq. 
(I5T1) and (1521) . 

The simulations of Figs. [T5] (inset), [15] and [16] have 
been made for i = lm and w c = 27t x 10 9 s -1 . These 
parameters are those of a realistic experiment on water 
and, such that the Brillouin precursor generated by the 
canonical input signal predominates on the main field 
which remains visible. The analytical results are, how¬ 
ever, general. They notably show that the validity do¬ 
mains of the low- and high-conductivity approximations 
are not entirely determined by the absolute value of the 
conductivity as considered in Ref [0 but strongly de¬ 
pend on the propagation length. For example, a conduc¬ 
tivity a ~ 3.3 x 10~ 2 S/m suffices to be perfectly inside 
the low conductivity domain (i? ss 51) when < = lm, 
whereas a conductivity as low as a m 1.72 x 10~ 3 S/m is 
required to attain the same value of R when £ = 100 m . 
Note additionally that the amplitude-reduction factor F 
is as large as 38 in the latter case but 2 in the former. 


VII. CONCLUSION 

Sommerfeld and Brillouin precursors originate from 
interrelated effects of group-velocity dispersion and of 
frequency-dependent absorption of the medium. The for¬ 
mer prevail in the formation of the Sommerfeld precur¬ 
sors jl], ii but both effects generally contribute to 
that of the Brillouin precursors 010. In Debye me¬ 
dia, however, as soon as the Brillouin precursor is dis¬ 
tinguishable from the main field, we have shown that 
the impulse response of the medium is perfectly fitted 
by the Gaussian obtained by neglecting the dispersion 
effects. By means of simple convolutions, we have then 
obtained explicit analytical expressions of the response of 
the medium to a wide class of input signals and extracted 
from them those of the precursor. Our main results are 
summarized below. 

For input signals with a single initial discontinuity, the 
Brillouin precursor is well approximated by a Gaussian 
(Figs. [51 and [Til) or a Gaussian derivative (Figs. [5]and[5]), 
depending on the parity of the discontinuity order. This 
approximation is excellent in the case of the canonical in¬ 
put signal uh {t) sin (u> c t) considered by Sommerfeld and 
Brillouin (Fig. [3]) and an exact expression of the precur¬ 
sor amplitude is obtained in this case [Eq. (fT51) |. When 
the input signal has two successive discontinuities sepa¬ 
rated by one period of the carrier frequency, the quasi¬ 
destructive interference of the corresponding precursors 
may generate a unique precursor with a second-Gaussian- 
derivative shape (Figs. 171 fH/li. and fTSl) . 

In the limit where the optical thickness of the medium 
at the carrier frequency is very large, the above- 
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mentioned shapes of precursor are exact and the precur¬ 
sors are simply proportional to the impulse response, to 
its first or to its second derivative, with amplitude scaling 
with the propagation distance £ as £~ 1 ^ 2 , £~ x and £ -3 / 2 
respectively. These remarkable asymptotic properties are 
not specific to particular input signals but are common 
to the precursors generated by all the input signals hav¬ 
ing the same integral properties. The same precursor can 
thus be generated by quite different input signals (Fig. 

m- 

An eventual static conductivity of the medium does 
not affect significantly the shape of the output signal 
and in particular of the precursor when the parameter 
R given Ea. (l53l) is large (low conductivity limit). The 
main effect of the conductivity is then an overall reduc¬ 
tion of the medium response [see Eq. (1481) and Fig. [H| . 
In the opposite case where R <Cl (high conductivity 
limit), the precursor disappears and the output signal is 
well approximated by the broad signal obtained by ne¬ 


glecting the polarization contribution, the amplitude of 
which scales as £~ 2 [see Eq. (15^1) and Fig. US] . 

The results obtained in our study cover various sit¬ 
uations. Only the Brillouin precursor generated by the 
canonical input signal has been actually evidenced in real 
experiments [27]. Our theoretical analysis is expected to 
stimulate complementary experimental works. 
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